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Abstract 



We show that two tight binding electrons that repel may form a bounded pair in two dimensions. 
—I The paired states form a band with energies that scale like the strength of the interaction 

Q^ potential. By applying an electric field we show that the dynamics of such states is that of a 

J, composite particle of charge 2e. The system still sustains Bloch-like states, so that if the two 

5/2 bands overlap single and paired states might coexist allowing for a bosonic fluid component that, 

if condensed, would decrease the resistance at low temperatures. The presence of two bands 

allows for new oscillations whose experimental detection would permit a direct measurement of 

the interaction potential strength. 
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O I. Introduction 
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Bound states are normally associated with basins of attraction, and it is always surprising to find 
them in the presence of a repulsive interaction. The most remarkable case is the pairing of electrons 



t" — in normal superconductors, mediated by lattice distortions. A new kind was recently observed for 



ultracold rubidium atoms in an optical lattice [1], where the binding arises from pure quantum 
interference. The pairs were shown to be stable, thus suggesting that they might form a superfluid 

1/^ phase if the composite is a boson. 

^^ The effect was originally predicted for a one dimensional string [2-5]. The eigenstates decay 

exponentially in the relative coordinate, while in the center of mass coordinate they are extended 
over the whole lattice. In higher dimensions one can argue that paired states may exist as well. 
To see this, consider two particles that repel moving in a two dimensional lattice. Together, the 
pair combine four spatial degrees of freedom, thus allowing that it be formally described as a single 
^ particle moving in a 4D lattice. A particle-particle interaction appears then as an interface potential 

in four dimensions, which decays to both sides of the hyperplane xi= X2, yi= y2, where (xi,yi), 
(x2,y2) are cartesian coordinates for the particles position in real 2D space. Surface states bounded 
to the interface are then expected even if the potential locally rises, as would occur under particle- 
particle repulsion. Since the interface represents matching coordinates for the two particles in 2D 
such surface states are paired states in the lattice. A similar argument applies in 3D. 

We find the 2D case particularly interesting since the onset of pairing could be relevant to high 
temperature superconductors where transport is related to conducting sheets, and to certain 2D 



experimental probes for which the physics is stih controversial [6-16]. It is believed that in the latter 
case localization due to disorder on the conducting sheet should produce a divergent resistivity as the 
temperature is lowered, but this behavior was not observed at all electron densities in metal-oxide- 
semiconductor field-effect transistors (MOSFET). In fact, a metal insulator transition was found 
in low density 2D samples as the density is slightly increased [8]. An early prediction invoked the 
strongly interacting nature of a low density electron fluid [9-10], while more recent theories study 
the effect of formation of a Wigner solid [13] or a superconducting phase [14], a strong spin-orbit 
interaction [15] and a few classical effects [16], none of which has received consensus for describing 
the physics at the root of the transition [17]. 

As we shall show, paired singlet states are predicted to exist in 2D for repelling particles in a 
lattice. In Section II we define the model and show that the density of states reveals the presence 
of two bands of different character, one corresponding to a Bloch particle in 4D and one describing 
a 2D surface state. In Section III we show that the dynamics of such surface state is that of a 
composite boson of charge 2e, our main finding, and in Section IV we discuss our results. 

II. The Model 

Consider for definiteness two interacting tight-binding electrons in 2D in a uniform external field. 
The asociated Hamilton operator for such system reads: 

-" ^ -^ / , / , \'^x,y+l,s'^3:,y,s + '^x,y,s'^x,y+l,s + '"'x+l,y,s'^x,y,s + '^x.y.s'^x+l.y.s) 



x,y s 

n 

x,y s x,y 



e^a— /|- X! X! ^l,y,sCx,y,s + UY^ c\y^^Cx,y,^c\y^^Cx,yX, (1) 



where cj. „ ^ and c^^y.s are the creation and annihilation operators for one electron located at the 
site [x^y] with spin s =tiii -^ the usual hopping energy parameter, n the sum of position vectors, 
e the charge of the electron, E and n the electric field magnitude and direction, respectively, and a 
the lattice parameter. The last term represents a two-body Hubbard contact interaction potential 
\J . We assume that the interaction is strongly screened and may be ignored for particle-particle 
separations beyond a lattice constant and discuss only the singlet state [3]. 

In the Wannier representation, the time dependent wave function for the pair in a singlet state 
may be expanded as 

l<^W)=X! X! •^^i'!'i;^2,!/2 Wki,2/i,si;a;2,y2,S2) (2) 

xi.yi X2,y2 

where the sum runs over all lattice sites, fxi.yi;x2,y2 (^) i^ the amplitude for an electron with spin si 
to be at (xi,yi), and another with spin S2 ^ si to be at (x2,2/2), while the ket \xi,yi, Si;x2,y2, S2) 
represents such state. The amplitude obeys the Schrodinger equation of motion 

*""77/a;i,j/i;a;2,J/2 ^ ^ ^Jxi,yi;X2,y2 < \EtjaA. + '-^ '^x-i,X20yi,y2) Jxi,yi:X2,y2 J \'^) 

where / is the sum of amplitudes of all nearest neighbors to the site {xi,yi^X2,y2) in the lattice, 
and X = {xi +yi+ X2 + y2)/\/2 in units of the lattice constant. We have assumed that the electric 
field is along the diagonal xi = j/i, X2 = j/2- 



By making the replacement 

the associated eigenvalue equation in the absence of an external field becomes 

Eg{u, v) = — 2A [cos k^a {g{u — 1, w) + g{u + 1, v)) + cos kya {g{u, w — 1) + g{u, v + 1))] 

+ Su,oSv.oUg{u,v) 



(4) 



where now u ^ X2 — xi,v — 1/2 —^ j/i.One can easily check that g{u, v) — Sufi^v,o is an eigenvalue of 
this equation at the lowest band edge kx — ~ ky, oi eigenvalue E — U. This solution represents 
extreme pairing, when one particle is precisely on top of the other. 

In the absence of interaction solutions of Eq. (4) are plane waves in a Bloch band, with eigenvalues 
E = — 4A(cosfca;a + coskya), where k — {kx,ky) is the center of mass wave number. The band 
remains when the repulsive interaction is turned on, while additional states appear at positive 
energy. This is seen in Fig. 1, where the density of states obtained by numerical diagonalization of 
the eigenvalue equation for several positive values of the parameter U is shown. Data corresponds 
to a finite lattice with N=49 plaquettes using periodic boundary conditions. The bell shaped curve 
represents the density of states of a single particle in 4D, the analog problem in this case, while a 
2D-like density of states profile representing surface states in 4D moves away as U grows. Note that 
the lower edge of this band is at E — U as expected. 




Figure 1: Density of states of two interacting electrons in a 2D lattice of 49 sites. For U — 0, it is 
that of a single tight-binding particle in 4D. When C/ 7^ a new subband emerges, moving to higher 
energies as U is increased. 



III. Evidence for bounded pairs 

To confirm the presence of two different kinds of states, we fiave calculated the average distance 
between particles 



d-i 



a;i,2;2 yi,y2 



(a;i-X2) +{yi-y2) |/;^i 



yi\x2,v2 I 



(5) 



on a finite 2D sample of N plaquetes using periodic boundary conditions. Here j labels the eigenstate 
considered. The results are shown in Fig. 2, where the quantity dj is shown as a function of sample 
size both for states in the B-band as well as states in the U-band. We see that the average separation 
of the particles is essentially null for states in the U-band for all cell sizes, as is expected for a bounded 
pair. For states in the Bloch band the separation is finite, with a value that grows linearly with N^ 
as expected for extended states in two dimensions. 
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Figure 2: Average mean distance of Bloch and U states plotted for different lattice sizes and a 
contact potential U = 50A. The mean separation in a Bloch state (B-band) grows linearly with the 
edge size, while for the U-band it stays near zero as expected for paired states. 



A more dramatic confirmation that repelling electrons in a lattice form bounded pairs is found 
when the electric field is turned on and the time evolution of a given initial state is followed by 
solving Eq. (3). It is well known that charged particles in a lattice experience Bloch oscillations 
in the presence of the field, with a frequency proportional to the particle charge[18]. To find out if 
such oscillations are present we studied the time evolution of the average position along the x axis 
of each particle 

i?ni)=«E E \f-^m■,-2,V2{tf X^, (6) 

xi,yi X2,y2 

where i — 1,2, is the particle index, and a similar equation for the coordinate y. The numerical work 
was done using a half implicit numerical method which is second-order accurate and unitary [19]. 
The positions were indeed found to oscillate, not simply performing Bloch-type oscillations but a 
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Figure 3: Frequency spectrum of the average position of each particle in a 2D lattice of 81 sites and 
an external field of strength F = lOA. (a) is in the absence of interactions and (b) is for a Hubbard 
interaction strength U — lOOA. The inset shows details near the Hubbard frequency. Frequencies 
are in units of A. 



more complex pattern depending upon initial conditions. Figure [3] shows the power spectrum for a 
2D lattice of iV = 81 sites and an electric field energy F — eEa = lOA, set rather large in order to 
avoid reflections from the edges in our reduced numerical sample. In the absence of interactions only 
the Bloch frequency w — F {h — 1) has any weight, as seen in Fig. 3(a). We next set the interaction 
strength at f/ = lOOA, a large value appropriate to have a U-band well separated from the extended 
states. For an initial condition with zero amplitude save for points far from the interaction region 
xi — X2, yi — y2, the situation is unchanged. But for an initial condition with unit amplitude 
in sites over the interaction hyperplane in the 4D space of coordinates (a^i, yi, 2^2: 2/2) = i''',s,r,s), 
r, s = 4, 5, 6, and either r or s = 5 when r ^ s, representing a portion of the hyperplane and a few 
nearby points, we see a Bloch oscillation with twice the Bloch frequency, with a power spectrum 
as shown in Fig. 3(b). We interpret this new frequency as arising from motion of a composite 
particle of charge 2e in an electric field, thus confirming that the pair behaves dynamically as a 
stable composite in the absence of dissipation. Additional interaction-induced oscillations (ININO) 
appear near U + F and U — F, shown in more detail in the inset. These latter oscillations, seen also 
in ID systems, are the result of pure correlated electron transport [20]. We note that our numerical 
results show that the features just described persist as the interaction strength enters the region 
J7 ^ A, only that hybridization of localized and extended states makes the power spectrum more 
complex. 



IV. Discussion 

In summary, we have shown that two interacting electrons moving in a lattice in 2D have singlet 
paired states grouped in a band, with the dynamics of a composite particle of charge 2e. Our 
results were obtained for a highly screened contact repulsive potential, yet studies in ID show that 
a Coulomb tail does not destroy the pairing [2] in the singlet state. Although the specific treatment 
refers to electrons in a lattice, the main results may be extended to other systems such as atoms in 
an optical lattice [1]. 

In order to understand the origin of the paired states we recall that ordinary Bloch states may 
be thought of as atomic states that hybridize with lattice neighbors, forming a tunneling network 
capable of sustaining extended states that are grouped in bands. Similarly, the paired eigenstates 
may be described as a high energy two-electron single well state that forms a band of extended 
states owing to pair tunneling throughout the lattice. The Hubbard model captures the escence of 
this picture, but details such as possible molecular two-electron extended states may arise if a finite 
range interaction is included. 

The energy of the paired states scales like the interaction strength U. If this latter quantity is 
of the order of the bandwidth, such states may lie partly or wholly within the band of extended 
states, allowing pairs to form in the ground state if there are enough electrons in the sample. The 
paired states equal in number the extended states, since they correspond to single particle surface 
states over a 2D planar interface in a 4D lattice. We thus speculate that if condensation of such 
pairs indeed occurs then a transition to low resistance should take place as the density is increased, 
as is actually observed in experiment [11]. 

Our model yields a pair whose dynamics is that of a stable composite of charge 2e. Such stability 
may be reduced by dissipation mechanisms such as phonon excitation. In a static periodic potential 
such as that generated by an optical lattice phonons are not present, and a pair of energy ^ U > 8X 
above the band center is expected to be stable since its decay into a band state is forbiden by energy 
conservation. Under an external field F, transitions between extended and localized states give 
rise to oscillations. If the associated frequencies to^ = U + F and a;_ = U — F were accessible to 
experiment, then the sum aj++ w_ = 2U would yield a direct measure of the interaction strength 
in the 2D sample. 
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